Relativistic Energy Levels 
of Para-Helium 



R. Grabeldinger, P. Schust, M. Mattes and M. Sorg 
II. Institut fur Theoretische Physik der Universitat Stuttgart 

Pfaffenwaldring 57 
D- 70550 Stuttgart 
Germany 

e-mail: sorg@theo2 . physik . uni-stuttgart . de 
Abstract 



The practical usefulness of Relativistic Schrodinger Theory (RST) is tested by calculating approx- 
imately the energy difference between the excited singlet state ls2s 1 So and the ground state Is 2 1 So 
of the helium-like ions with arbitrary charge number z cx (2 < z cx < 100). The results are compared 
to the corresponding predictions of other theoretical approaches in the literature and to the experi- 
mental data. Since the exact solutions of the RST energy eigenvalue problem are unknown, one has 
to resort to approximative methods. However the crudest approximation ("spherically symmetric 
approximation") yields relatively accurate results so that it seems worth while to develop more 
powerful approximation techniques 



I. INTRODUCTION 

The well - known wave - particle duality [lj in quantum theory Q enables one to interprete 
some of the quantum phenomena in terms of the particle picture ( with the associated statis- 
tical interpretation), whereas other phenomena of elementary matter are more conveniently 
understood in fluid - dynamic terms due to the wave picture, e.g. the Bose - Einstein conden- 
sates 0,0]. Especially interesting are the fields of "intersection" where both pictures can be 
simultaneously applied, and then must be expected to predict the same numerical results for 
the outcome of the corresponding experiments. An example for this type is non- relativistic 
atomic physics where it can be shown that the fluid - dynamic approach (i.e. the density 
functional formalism 0,0]) is exactly equivalent to conventional quantum mechanics, being 
based on the use of wave functions in place of density distributions. 

The situation is somewhat different in relativistic quantum mechanics where neither a 
well -working particle theory nor its fluid - dynamic counterpart for N- particle systems does 
exist ( see the many critical comments on the Bethe- Salpeter equations, for instance in ref.s 
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The recently established Relativistic Schrodinger Theory ( RST ) jtj Q is intended to fill 
just this theoretical gap where both basic concepts ( i.e. wave functions and densities ) be- 
come united into one picture. From the mathematical point of view, this new theory is 
related rather to the fluid - dynamic than to the probabilistic view upon quantum theory. 
The reason for this is that RST takes the Whitney sum of the single - particle bundles as the 
many - particle bundle, in contrast to conventional quantum mechanics which is based upon 
the tensor product of the single - particle Hilbert spaces. Subsequently, we will present 
an example where the fluid - dynamic approach (i.e. RST) is able to predict frequencies 
of atomic spectral lines with acceptable precision, albeit somewhat more inaccurate than 
conventional quantum mechanics does. It is important to place emphasis on the fact that 
we have to compare here our approximative RST results to those of other approximative 
methods, which have to be applied in the conventional domain because a consistent and 
complete conventional theory for relativistic N - particle systems does not exist. More con- 
cretely, we will apply RST in order to compute the energy difference between the excited 
singlet state Is 2s x 5o and the ground state Is 2 of the helium -like ions with arbitrary 
nuclear charge number z ex ( 2 < z cx < 100). Furthermore, we will compare our results to 
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the analogous predictions obtained by the 1 jZ - expansion method which are close to 
the analogous predictions of the all -order technique in relativistic many- body perturbation 
theory ( MBPT ) [12j. Since the latter two approximation techniques rely on the conventional 
tensor product formalism, whereas RST utilises the Whitney sum, the predictive po- 
tentiality of both formalisms can be directly opposed. And it is worth noting that RST's 
fluid - dynamic approach delivers predictions relatively close to the observational data; see 
tables I, II and figures 1,2. However, it must be stressed that the present comparison of 
numerical predictions for spectroscopic data between RST and the conventional approach is 
actually a competition between approximative methods. Therefore, the fact that the present 
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RST results are somewhat more inaccurate than their conventional counterparts 
not hint on an intrinsic deficiency of RST itself but merely points to the necessity of devel- 
oping better approximation techniques within the framework of RST. Such an improvement 
may be conceived in a two - fold way, namely by 

(i) invention of more adequate trial functions for extremalising the RST action integral 

(ii) solving the corresponding two - dimensional mass eigenvalue problem in the variables 
r, $ (in place of resorting to the spherically symmetric approximation). 

The results are elaborated through the following arrangement: 

For the convenience of the reader unfamiliar with RST, some of this new theory's funda- 
mentals are collected into a brief survey, such as the RST dynamics in section [H] and the 
conservation laws following from the RST dynamics in section IIIII But of no lesser impor- 
tance is the RST kinematics, such as the fibre metric K a p for the Lie algebra bundle. Indeed, 
this object establishes an important link between the RST currents j afl and the Maxwell 
currents j 01 ^ emerging in the gauge field equations; see sections HVI andlVl The existence of 
an action principle ( section IVII ) is useful not only for the logical coherence of the theory 
and for the deduction of the conservation laws by means of Noether's theorem, but also for 
establishing approximation methods based upon variational techniques. The point here is 
that the RST eigenvalue problem for computing energy levels of bound N - particle systems 
is sufficiently complicate in order to prevent elaboration of analytic solutions. Therefore, 
we resort here to the spherically symmetric approximation. This method consists in adopt- 
ing spherically symmetric ansatz functions for the wave amplitudes and substituting this 
"isotropic" ansatz into the action principle in order to deduce the corresponding 50(3)- 



symmetric eigenvalue equations in section IVIII and lVIIII There exists an interesting point 
with this isotropic ansatz; namely, it is exact (!) for the states ns 2 ^o, which share their 
symmetry with the ground- state Is 2 1 Sq. The exact wave amplitudes of the more general 
solutions rtis n 2 s x Sq with different principal quantum numbers (fi] ^ n 2 ) are not spher- 
ically symmetric but depend on both spherical polar coordinates r and However, the 
isotropic ansatz with its exclusive dependence of the wave functions on the radial variable 
r is tentatively considered as an acceptable approximation for those non - isotropic states 
(m ^ n 2 ). 

The exactly isotropic states ns 2 were already considered in great detail in some of 
the preceding papers and are therefore presented here only briefly in section 11X1 mainly in 
order to introduce the RST concept of self- interactions which is based on the fibre metric 
K a p. The latter object is fixed up to one free parameter u (the self '- interaction param- 
eter) whose determination has to be performed by adjusting the RST predictions to the 
corresponding experimental value for one selected physical situation (see the selection of the 
bismuth ionisation energy in [lfll l"!^). 

The spherically symmetric ( and therefore approximative ) eigenvalue equations for the 
states riiS n 2 s So with n\ 7^ n 2 are presented in section [Xj followed by the normalisation 
prescription for the corresponding solutions in section IXII The occurrence of the self- 
interactions is much more manifest in the general case (ii! / ri 2 ) than it is in the exactly 
isotropic subset ns 2 So of section IIX1 The latter subset is easily identified within the more 
general set n\s n 2 s by letting coincide the principal quantum numbers (ni = n 2 ), mass 
eigenvalues (Mi = M 2 ), gauge potentials, etc. In this sense, the exact isotropic states 
ns 2 Sq are kinematically included in the larger but approximative set n\s n 2 s Sq. 

Finally, for the theory to be able to produce definite predictions, it needs an energy 
functional (Et, say) whose value upon the numerical solutions of the exact [n\ = n 2 ) or 
approximative ( n% 7^ n 2 ) eigenvalue system yields the desired atomic energy levels. This 
RST energy functional is the subject of section IXII| and the corresponding numerical re- 
sults are displayed in section IXIII| see tables I, II and figures 1,2. The corresponding 
experimental values are taken from ref. [l^ |. 
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II. RST DYNAMICS 



The RST dynamics subdivides into the matter dynamics, Hamiltonian dynamics, and 
gauge field dynamics. The matter dynamics consists of the Relativistic Schrodinger 
Equation ( RSE )for the wave function \I/ 

ihcv^ = , (i) 

or, resp., the Relativistic von Neumann Equation (RNE) 

v,i = -jr [in, - n,i] , (2) 

if matter is to be described in terms of an intensity matrix I rather than in terms of a wave 
function Clearly, a pure state \I> can always be conceived as a special type of mixture X, 
namely, as the tensor product of \I/ and its Hermitian conjugate 

l^^(g)^. (3) 

The Hamiltonian 7i M is itself a dynamical object of the theory; and therefore it is to be 
determined by its own field equations, i.e. the integrability condition 



and the conservation equation 
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- j-WH, = -the | ^— j + YTT^ ) . ( 5 ) 

The meaning of the integrability condition (0]) is to ensure the bundle identities, such as 

\VpV v - V V V^\ * = Tpft (6) 

[V ll V v -V u V IM ]Z= , (7) 

whereas the conservation equation (|3J) guarantees the existence of conservation laws, for 
instance for the total current 

V% = . (8) 
Such a conservation law is indispensable for defining the particle number N through 

N = [ j^dS" , (9) 



where the choice of the hypersurface (S) is arbitrary, just on behalf of the local law (EI)! 



III. CONSERVATION LAWS 



The matter densities, such as the total current or the energy- momentum density 
^Xr^, are quite generally defined by means of the intensity matrix I and a corresponding 
operator (e.g. velocity operator or energy - momentum operator T^ v ) by the following 
trace prescription: 

^ = tr(x.r M ) (10) 

^% v = ti{X-% p ) . (11) 

For the special case of N identical Dirac particles, the total velocity operator V M is taken as 
the iV-fold direct sum of the ordinary Dirac matrices 7^ 

r M = 7m © 7m © • • • © 7m , (12) 

V v ' 

at -fold sum 

and the energy - momentum operator is assumed to be of the following form 

V = 4 {T»Hv + HuTf, + YyU^ + H^Yy] . (13) 

If the considered N - particle system is acted upon by an external field ^F^, the energy - 
momentum tensor of matter develops a non- trivial source which equals the Lorentz 

force density ( L >/„ ^ Q 

= W/„ . (14) 

Since this force density is composed of internal and external contributions, it does not 
vanish for vanishing external field strength ^F^ v . However, if one includes also the energy - 
momentum density of the gauge fields ^■T^ and the external interaction density in 
order to obtain the total energy - momentum density ^^T^ as 

{Ty T^ = (D K» + {G K + , (15) 

then this total density has vanishing source whenever no external force field ^'F^ due to 
an external source ( - ex - ) j A1 is present (i.e ^ ex F^ = 0, ^j^, = 0): 

= hcF,Jrtf ^ . (16) 

In contrast to this energy -momentum conservation, which does hold exclusively for closed 
systems only ( = ), the charge conservation law ((HI) is valid in any case, whether the 



system is energetically closed or not. This can easily be demonstrated for the pure states 
(01) when the following alternative version of the conservation equation (01) is used js, 10] 



T»H^ = Me 2 , (17) 

which directly relates the Hamiltonian TC^ to the mass operator M.. Indeed, multiplying 
both sides of the RSE (U) by the total velocity operator T M from the left and observing 
( IT7I ) yields just the Dirac equation for the iV- particle wave function \1/ 

ihVV^ = Mem . (18) 

Furthermore, the total current j M ( fTUl ) reads for the pure states $ (01) 

j, = , (19) 

so that the divergence of vanishes just on account of the Dirac equation (EH) (and its 
Hermitian conjugate) where the covariant constancy of the velocity operator T M has also 
been presumed 

V^T U = (20) 



IV. CURRENTS 



In order to close the RST dynamics, one has to supply some field equation for the bundle 
connection A^. This object enters the RST dynamics via the covariant derivative £> M of the 
field objects, e.g. for the wave function \1/ in the RSE (^) 

V^ = d^ + A^ , (21) 

or similarly, for the Hamiltonian Ti^ in the integrability condition (0]) 

V^Hv = Vfj,H u + [Ap, H v ] ■ (22) 

The choice for such a field equation for the "gauge potential" A^ is not quite arbitrary 
because the desired equation has to be compatible with the already fixed RST dynamics. 
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In order to exemplify this important point for a special RST object, consider the gauge 
currents j afl , (a — 1, ... , N 2 ) 

jan = tr (J • Van) > ( 23 ) 

with the gauge velocity operators v afJ- being defined in terms of the gauge algebra generators 
r a through 

v ail = l -{r a ,T ll } . (24) 

It can be shown that the RST dynamics implies the following source equations for these 
gauge currents j afM (1231): 

D»j*n = , (25) 
with the covariant derivative being defined through 

DfJaa, = V^jav - UJ afl jp v . (26) 

Here, the connection one -form = \uj^ a n\ takes its values in the adjoint representation 
of the original connection 

An = A V* , (27) 

i.e. one has 

with the U(N) structure constants C /3 7Q being defined as usual: 

[r Q ,r /3 ] = C\ /3 r 7 . (29) 
The source relation (EH) may also be recast into the more compact form 

V» J„ = , (30) 
when the current operator is defined as follows: 

Jn = y>« ■ (31) 

Here, we have used the contravariant current components ("Maxwell currents" ) which 
are linked to their covariant counterparts j afl ( " RST currents" ) by means of some fibre 
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metric K a p {compatibility tensor) 

j\ = K^jfr (32) 
Jctfi = K aP f^ (33) 
(K af3 K^ = 5\) , 

which itself has to be covariantly constant: 

D^K a/3 = . (34) 

Indeed, this covariant constancy of the fibre metric K a p is necessary in order to ensure the 
source relation ( 1251 ) also in contravariant form 

D^f, = , (35) 

which, in turn, acts as the desired integrability condition for the gauge field equations. Thus, 
the RST matter dynamics with its "conservation law" (EH) becomes compatible with the 
gauge field dynamics requiring the "conservation law" (13*51) . 

V. GAUGE FIELD EQUATIONS 

After all these arrangements, it is now almost self- suggesting to select a field equation 
for the bundle connection A^. For this purpose, first consider its curvature T^ v 

•Fp V = V^v4.j/ Vj/vA^ + [vA^j, J\.u\ , (36) 

and then postulate that this object obeys the ( non- Abelian ) Maxwell equations 

= -*ma s J v (37) 

( as - k) ■ 

In components, this gauge field equation reads 

= 4tt« s j^ , (38) 

provided the curvature J 7 ^ is decomposed into its components F a ^ u ( "field strengths" ) as 
usual; cf. ( 1271 ) 

•F l uv F fn/7~a • (39) 



Indeed, now one is easily convinced that the Maxwell equations ( 1371 ) are compatible with 
the preceding RST dynamics, and therefore also with their implication (1301) . Namely, it is 
just the generally valid bundle identity 

V^T^ = , (40) 

that has to be applied to the Maxwell equations (EH) in order to see that the continuity 
equation (1301) is actually satisfied whenever the Maxwell equations are obeyed. 



VI. ACTION PRINCIPLE 



The logical consistency of the coupled RST - Maxwell dynamics is supported by the fact 
that it can be deduced from an action principle 

(Wrst = , (41) 

where the RST action integral is given as usual in terms of a Lagrangian Lrst through 

W RST = J d 4 xL RST [^,^] . (42) 

The preceding conservation laws may then be deduced from this action principle via the 
Noether theorems. According to the subdivision of the RST fields into a matter part ^ 
and a gauge field part A^, the total Lagrangian Lrst will be also composed of a matter 
contribution and a gauge field part Lq, together with a contribution L ex due to the 
external fields (e.g. Coulomb potential of a nucleus): 



7 RST 



L D + L G + L ex . (43) 



Here, the matter part L-^ reads in terms of the wave function \1/ 

ihc 

L D = — [ tyr^V^ - {p$) M ] - Mc 2 ^ , (44) 

the internal gauge field contribution Lq is given in terms of the gauge potentials A a ^ as 

ftp 

l " K 07 F^ U F^ (45) 



lQrra s 
he 

16^' V ^' ' 
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and finally the external contribution L ex contains the external objects ^'A^, ^'F^, 
in the following combination: 

L cx = he ( K4 M • f - ■ A" - J—WF ia ,F> a '\ (46) 
\A-n ^ A 1 p + A 2 p , F^y = V fjiAn — V V A^ 



With this arrangement, the extremalisation of the action integral (021) with respect to the 
wave function \1/ yields the RST-Dirac equation (El); and similarly, the variation with 
respect to the gauge potential lets emerge the RST- Maxwell equation (1371). 

Clearly, the scientific value of a new theory has to become evident by its successful appli- 
cation to practical problems. For this purpose, we choose here the energy difference A£ , 1 _ 2 
of the para -helium states ls2s 1 Sq and Is 2 The reason for this choice is that, on the 
one hand, for this quantity there do exist sufficient experimental data Jl4j, an d on the other 
hand, the corresponding predictions of other theoretical approaches Il2l | are also avail- 
able. Thus, a comprehensive comparison of the predictive power of the various theoretical 
frameworks becomes feasible, both with inclusion or neglection of self- interactions. 

VII. STATIONARY BOUND STATES 

It should be a matter of course that the desired energy difference A£ , 1 _ 2 emerges also 
in RST by solving as usual the associated energy eigenvalue problem. The corresponding 
RST eigenvalue equations are to be deduced from the general RST dynamics by means of a 
stationary ansatz for the two -particle wave function \1/ 



*(*) = • (47) 




The ansatz for both Dirac spinors ip a , ( a — 1, 2 ) is as usual 

M a c 2 " 



i> a {r,t) = exp 



h 



-t 



^a(r) , (48) 



with mass eigenvalues M a to be determined from the stationary form of the two - particle 
Dirac equation (ITH1). For this purpose, the time - independent parts ip a (r) of the Dirac 
spinors (SHI) are further decomposed into two - component Pauli spinors ( a ty±(r) according 
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to 



1 1 



Xr) = R + (r)-(2> 2 (49) 

.(7=) = -i\R_(r) • C*i* ( 50 ) 

.(0 = £ + (r) • C*if * ( 51 ) 

.(r) = -i5_(r) -CH^ • ( 52 ) 



Here, the usual spinor basis C, 3 \ m |15( has been used, obeying the angular- momentum rela- 
tions 

Ht = j u + 1) »vr (53) 

J 2 C'| m = mKT (54) 
L 2 C^ m = Z (Z + 1) n 2 C^ m (55) 

s 2 c°r = s( s +i) wr » (56) 

with the addition theorem for angular momentum 

3 =l±s. (57) 

For electron spin s — ~ this yields the configurations with j = ~; Z = 0, 1; m = ±~ 
when restricted to singlet states ^o- It turns out that the presumed spherical symmetry 
("isotropic ansatz" ) of the wave amplitudes R±(r), S±(r) is an exact symmetry for the 
states ns 2 1 S$. But for the general set n\s ri2S 1 So, this symmetry is adopted to be a useful 
approximation. 

In a similar way, the gauge fields A a M have to be specialised now to the stationary 
situation. 

VIII. GAUGE POTENTIALS 

First, the external potential ( ex l4 M = ^ cx ^o(r*); — A ex (f) J is assumed to be due to a fixed 
point -like nucleus of charge number z ex . Therefore, we put 



K4 (r) = , (58) 



with no external magnetic field being present (i.e. A ex (r) = 0). Next, the electromagnetic 
potentials A a ^, a = 1, 2, cf. (123), are also assumed to be 50(3) or 5*0(2) symmetric, resp., 
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and static: 

A%^{^A {r);-A a (r)} , (59) 

where both electrostatic potentials ^ a 'Ao(r) are additionally assumed to be identical for the 
special subset of isotropic states ns 2 ^So: 

«4)(r) = ( 2 U (r) = ^Uo(r) . (60) 

Since the spins point in opposite directions for para- configurations singlet states), the 
magnetostatic vector potentials are assumed to differ in sign for this special situation: 

A 1 {r) = -A 2 {r)=A p {r). (61) 

It is self - suggesting to evoke such a structure also for the electromagnetic currents j afl 
(a = 1,2), since the Maxwell equations (EH1) link them to the electromagnetic potentials, 
i.e. we put 

= ^27^2 = k 2/I = { (2) k {r); -k 2 {r) } (62) 
J2 M = ^iTmV'i =f fa? = { (1) k (r); -kt(r) \ . (63) 



Again, for the special subset of states ns 2 ^So the electrostatic charge densities ( a -fc ( r ) do 
coincide 

(%(r) = {2) k (r) = {p) k (r) , (64) 
and the three - currents are antiparallel 

fci(r) = -k 2 (r) = k p (r) . (65) 

The simplest form of the three - current k p (r) and vector potential A p (r), to be used in the 
following for the singlet states ns 2 ^So, reads in terms of spherical polar coordinates ( r, </? ) 

k p (r) = k p (r) sin?? • e v (66) 
A p {r) = rA p {r) sin •& ■ . (67) 
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Clearly, for the more general states n\s n 2 s ^o, the currents and potentials will have 
different strengths ( k\ ^ k 2 , A\ ^ A 2 ): 

kiif) = ki(r) sin^e^ (68) 

k 2 {r) = k 2 (r) sin^e^ (69) 

A(f) = rAi{r) sintfe^ (70) 

A 2 {r) = rA 2 (r) sin-tfe^ , (71) 

but nevertheless, the scalars ( a 'ko, k a , ^A , A a are assumed to be still of the spherically 
symmetric form. The latter assumption spoils the exactness of the solutions for rii ^ n 2 
which can hold exclusively for the special subset n\ = n 2 { the exact but non - isotropic 
states riis n 2 s x Sq require angular dependent ansatz functions and potentials: R±(r, 
S + (r, $), etc., see ref. (l6| ). 

The exact 50(3) symmetry of the ns 21 So configurations also suggests the following shape 
for the magnetic exchange potential B{r) and associated current h(f) 

B(r) = irB(r)W p (r) (72) 
h(r) =ih(r)W;(f) , (73) 

where the complex - valued three -vector field W p {f) is given by 

W p (r) = -e iv (Ze + ico8&-e lp ) . (74) 

Here, the interesting point with the exchange four- potential B p = |.Bq, — Bj is that the 
ground -state symmetry (i.e. "isotropy" ) of the ns 2 ^So states demands the following iden- 
tity of the magnetostatic potential A p {r) ( IBTI ) and magnetic exchange potential B(r) (172*1): 

B(r) = A p (r) . (75) 

For the exact non - isotropic states nis n 2 s x Sq , the ansatz ( 1721 ) is retained, but now 
the assumption of spherical symmetry (i.e. putting B(r,d) =>- B(r) )spoils the exactness of 
the solution. Moreover, it is a general peculiarity of all singlet states that the electric 
exchange potential Bq(t) is zero, which is a consequence of the vanishing of the exchange 
density ho as the time- component of the exchange current h^ 

K ="= ^i7^2 = ~e 2u j\ = e 2u f% , (76) 
14 



with 

ho(r) = B (r) = . (77) 

This peculiarity of the para -helium states implies that there is no exchange energy of 
the "electric" type which leaves the occurrence of the exchange energy to its "magnetic" 
counterpart being based upon the space parts B, h of B^ and h^. Since these space parts 
B, h are usually much smaller than the time components Bo, ho, it may appear that RST 
( and the Hartree- Fock approach as its non- relativistic limit ) produces too inaccurate pre- 
dictions of the para -helium level system; see the discussion of this point in ref. Q|. The 
reason for this common deficiency of both the Hartree- Fock approach and RST is the fol- 
lowing: this refers to the approximative assumption that the two -particle wave function 
could be approximated by only two (albeit appropriately chosen) one -particle wave func- 
tions (see the preceding RST ansatz being specified by equations (El) - (EH))- Within the 
HF approach, the deficiency is superseded by resorting to the multiconfiguration Dirac- 
Fock method (MCDF) whose counterpart for curing the analogous RST deficiency will be 
presented in a separate paper. 



IX. MASS EIGENVALUE EQUATIONS FOR ni = n 2 

With the stationary form of all the RST fields being fixed, the next step is to substitute 
these expressions into the general two -particle Dirac equation (EH)- But only for the ns 2 ^So 
states, where both electrons share the same principle quantum number (ni = n 2 =^n), the 
present "isotropic" ( i.e. spherically symmetric ) ansatz can be an exact solution. For this 
special class of solutions, the following identifications are self-evident: 



Mi 


M 2 : M" 




R+{r) 


= S + (r) 




i2_(r) 


= S-(r) 




{1) k (r) 


= {2) Ur) = 


(p) k (r 


fci(r) 


= -k 2 (r) = 


fcp(r) 


Ai(r) 


= -A 2 (r) EE 


B(r) 



(78) 
(79) 
(80) 
(81) 
(82) 
(83) 
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The corresponding eigenvalue system for the wave amplitudes R±(r) is exact and reads 

d H (r) M 4- M" 

r ' + [ {ex Uo(r) + (p U (r)l • R(r) + 2rB(r) ■ R+(r) = - V c ■ RJr) 
dr L J n 

(84) 

d R (r) 2 M — M" 

d ~ l j + -#-(r) - [ (cx U (r) + (p U (r)] • i? + (r) - 2r£(r) • R_{r) = c • R+(r) . 

(85) 

Because the wave amplitudes R±(r) couple to the gauge fields ^A (r) and B{r), the eigen- 
value system has to be closed by supplementing the Poisson equations for these gauge fields, 
which are of course to be deduced from the Maxwell equations (1371): 

(tT + ~f) {P) Mr) = 4W p) Mr) (86) 
\ ar z r dr J 



^S + ~) B(t) + 6BHr) (l - §f"fl(r)) = 4 TOs e"*^ . (87) 

Here, the gauge fields ^A^r) and B{r) couple back to the wave amplitudes via the charge 
and current densities ^k (r) (EH) and k p (r) (ESI) which read in terms of the wave amplitudes 
fl±(r) 

% w = m±^i (88) 

47T 

Mr) = m^mi) . (89) 

A striking feature of the current mass eigenvalue system (EH) - ( 1871 ) is the occurrence of 
the self- interaction parameter u in the last equation (1871) . Obviously, this is the only place 
where the self- interaction effect enters the eigenvalue system. It can be shown that this 
peculiarity is due to the fact that for the ns 2 1 S states both electrostatic charge densities 
coincide; see equation (EH)- More generally spoken, the self- interactions emerge in RST in 
connection with the covariantly constant fibre metric K a p (E2) - (EH)- This objects converts 
RST currents j afl to Maxwell currents j a which, in turn, enter the Maxwell equations (1551). 
Thus, the fibre metric K a p essentially determines the coupling of the field strengths F a ^ v 
to the RST currents j afl and thereby acquires the meaning of a coupling matrix. Its general 
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{K a0 } 



(90) 



shape for the currently considered two -particle system looks as follows [loj | : 

/ e u ■ sinh u —e u ■ cosh u \ 
— e" • cosh u e u ■ sinh -u 
-c- 

V o -e 2 " o y 

Evidently, this fibre metric K a p owns one degree of freedom which plays the part of a 
"renormalisation parameter" , because it continously changes the coupling strengths of gauge 
fields and currents. This is the reason why it emerges in connection with the original 
coupling constant «s = €- on the right-hand side of equation (1571). It is not subject to 
any constraints, and therefore, it can be used to let the RST predictions coincide with 
the observational data. The point here is that this matching of theoretical predictions and 
experimental data can be accomplished by fixing one value for u for all the considered charge 
numbers z ex , see ref Observe also, that in the Poisson equation ( 1551 ) the electrostatic 
coupling constant as is not renormalised by the self- interactions, in contrast to the magnetic 
case of equation (1571): as =>■ ase~ 2u . As a consequence of this peculiarity of the states 
ns 2 l So, the magnetic interactions are completely switched off for u — > oo, leaving the field 
configuration to be of purely electrostatic type. This field configuration is already known 
from ref. [3] as the electrostatic approximation. 



X. MASS EIGENVALUE EQUATIONS FOR ni / n 2 

It is necessary to stress that the eigenvalue system ( EH ) - ( 1571 ) for the ns 2 1 Sq states is 
exact and therefore can be deduced either from the RST variational principle ( 1411 ) by means 
of the isotropic ansatz for the wave functions ip a and the gauge fields A a ^ or this system 
can be alternatively deduced by directly substituting the isotropic ansatz into the original 
field equations for matter fH5T ) and gauge fields (E7I). For a detailed study of the ns 2 x Sq 
state see refs. j^l [3, 16, Generalising now the situation to the states n\s n 2 s x So 

with different principal quantum numbers n\ and n 2 , the isotropic ansatz can no longer be 
substituted into the RST field equations, because they do not admit those isotropic, well 
localised solutions for n\ ^ n 2 (albeit solutions, diverging at spatial infinity, may exist). 
Nevertheless, the isotropic ansatz does not become useless for this more general situation, 
but it is merely devaluated to an approximative solution, whose quality may serve to estimate 
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the influence of the neglected anisotropy. The corresponding approximate field equations 
are obtained by substituting the isotropic ansatz into the action principle ( 1411 ) and carrying 
out the variational procedure. This yields the following mass eigenvalue equations for the 
wave amplitudes R± and S± 



dR ± + [(ex)^ + (2 Uo] . R _ _ 2 {A2R+ _ 2BS+) = _M 1 +A£ c . R 



dr 



h 



3 

<! + - • R-— [< CX U + ( 2 U ] • i? + + ^r (A a JL. - 2BS-) = Ml^L c . R 



dr r 



3' v ~~'" 7 ft 



dr 



3 

<lSL + - • S. - [H4 + (%] • 5 + - + 25i?_) = ^L^ c • 5. 



(91) 
(92) 
(93) 
(94) 



dr ' r l " "j 1 3 \ ^ 

In a similar way, the Poisson equations for the electrostatic potentials ^Ao(r) are found 
to be of the following form 



£ + is) (ll4 » w + r B » 



— - ( (1 U„(r) - < 2 U„(r)) 



4vra s • e~ 2u {K s (u) ■ ^k {r) - K p (u) ■ ( %(r)} (95) 



47ra s • e- 2u • <%(r) - tf p (u) • Mr)} . (96) 



Here, the exchange length parameter czm is defined as 

h 



iM — 



(Mi - M 2 ) c 
the pair coupling constant K p {u) is given by 

e 2u + 1 

KM = 

and similarly, the self - coupling constant K s (u) by 

2u _ i 
KJu) = — — . 



(97) 



(98) 



(99) 



Concerning the magnetostatic potentials A a (r) (1701) - (1711) and the exchange potential 
B(r) (1721). the analogous variational procedure yields the following ( generalised ) Poisson 
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equations 



\ ar z r ar J 

= 47T 
= 47T 



l--r 2 (Ax(r)-A 2 (r)) 



-2u 



{K s {u) ■ k 2 {r) - K p {u) ■ h(r)} 



l--r 2 (Ai{r) - A 2 (r)) 



-2k 



{K s {u) ■ h{r) - K p {u) ■ k 2 {r)} 



(100) 



101} 



/^d 2 _ + 4_d 
ydr 2 r dr 



B(r) + B(r) 



(Xm 



- (WA,(r) - (2 U (r)) 



n 2 



+ 



+3(Ai(r) - A 2 (r)) -r 2 



2 J B 2 (r) + -(A 1 (r)-A 2 (r)) s 



47rase 



-2u 



/i(r) 



with the current strengths k a (r) (1681) - (Ifi9"l) and /i(r) (ESI) being defined through 

R + (r) ■ R(r) 



k 2 {r) 
h(r) 



2tt 

S+(r)-S-(r) 
R + {r) ■ S-(r) + i?_(r) • S+(r) 

47T 



(102) 

(103) 
(104) 
(105) 



There is a nice consistency check for the more general situation with different quantum 
numbers n\ ^ n 2 : the previous eigenvalue system (EH) - (EH) for the special subset of states 



ns 



2 1 



So, being defined through the identifications (1781) - (1831). must be recoverable from the 
present system (1911) - ( fl02l ). Indeed, some simple arguments show that this requirement 
is satisfied. In this sense, the more general system for n\ ^ n 2 appears to be the unique 
isotropic ( and therefore approximative! ) generalisation of the special subcase with n% — n 2 . 
In view of this intimate relationship between the general and special cases, one may expect 
the general case yielding sufficiently realistic predictions for the ri\s n 2 s 1 Sq energy levels 
despite its approximative character; see the numerical results below for verification of this 
supposition. 

XI. NORMALISATION CONDITIONS 



The spectrum of desired solutions R±(r), S±(r) of the mass eigenvalue equations (EH) - 
) is unique only if a normalisation condition upon these solutions is imposed. Such a 
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condition is closely related to the asymptotic form of the electrostatic one - particle potentials 
( a l4o(r) (a = 1,2). Indeed, one expects that each of the two potentials ^Ao(r) is of the 
Coulomb form in the asymptotic region ( r — > oo ) 

K4 (r) =► , (106) 
r 

because each of the two electrons carries one negative charge unit. 

In order to see this asymptotic link between charge and potential more clearly, explicitely 
write down the first two equations of the RST - Maxwell system ( 1381 ) ( i.e. for (3 = 1,2): 

VFV = Anas {j\ - ^— [B^G% - B^G, V ] } (107) 

VF\ U = Anas |j 2 , + [B»G% - B^G^] } . (108) 

Here the two - particle connection A^ ( 1771 ) and its curvature (ESI) ar e decomposed with 
respect to the U(2) generators {r a } = {ti, r 2 ; x? x} as follows 

A„ = -* (CX U M • i + a\t x + a\t 2 + b^x - b;x (109) 

= -* {cx) i^ • 1 + F Vi + F\ v r 2 + G^ X - G;„x • (HO) 
Thus, defining the entanglement vector G M by 

G, = ^ [B v G* Vil - B*»G VIX ] , (111) 
both Maxwell equations ( 11071 ) and ( 11081 ) appear as 

= -47ra s (1) ^ (112) 
WF\ U = -4vra s (2 i , (113) 

provided the effective currents ^ a \ v (a = 1,2) are defined in terms of the entanglement 
vector G^ by 

(% = G„-j\ (114) 

{ X = ~G v -j\. (115) 

The motivation for introducing such an arrangement is that these effective currents must 
have vanishing source 

V M( % = 0, (116) 
(a = 1,2) 
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in order that the Maxwell equations ( 11121 ) - ( 11131 ) be consistent. As a consequence, the 
desired normalisation conditions can be imposed upon the wave amplitudes in the following 
form: 

J (%dS^ = 1 , (117) 

(5) 

where the choice of the hypersurface (S) in space -time is arbitrary just on account of 
the source equations ( 11161 ). The validity of these important source equations can easily 
be checked by straightforward computation: First, write down the first two non-Abelian 
source equations ( 1331 ) (i.e. for a = 1,2) 

VV, = i[BV M + 5"j 4 J (118) 

and then form the divergence of the entanglement vector ( II 111 ) under use of the definition 
of the curvature components F a ^ ( 13131 ) 

F\ u = V,A% - V v A a „ + C a Pl A\A\ , (120) 

which finally yields 

= i [B*»j\ + B"j\] . (121) 

Thus, combining the last three equations ( 11181 ). ( 11191 ) and ( 11211 ) actually yields the desired 
source equations ( 11161 ) for the effective currents ^ a \^ ( 11141 ) - ( 11151 ). 

Clearly, for the stationary bound states considered here, the hypersurface (S) of the 
normalisation integral ( 11171 ) will be chosen as a time -slice (t = const) which converts the 
general prescription ( 11171 ) to an ordinary integral over three -space 

J d 3 r (a) l (r) = 1 (122) 

(a = 1,2) . 

Furthermore, according to the general relations (13*2*1) . the Maxwellian currents j a v may be 
expressed in terms of the RST currents j afl ( "6*21 ) - ( 1631 ) such that the effective currents ^ a \^ 
( 11141 ) - ( 11151 ) finally reappear as 

{ \ = G U + e~ 2u {sinh u ■ k 2u + cosh u ■ k lu } (123) 
(\ = —G v + e~ 2u {cosh u ■ k 2v + sinh u ■ k lu } . (124) 
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Thus, denning the normalisation parameters z a (a = 1,2) as 

z a = I d 3 r^k (r) , (125) 



one finds that the normalisation conditions ( 11221 ) can be transcribed to the following con- 
straints for these normalisation parameters z a : 

h = ~ (1 - 9*) ■ K p {u) - (1 + g.) ■ K s {u) (126) 
z 2 = - (1 + g.) ■ K p {u) - (1 - g.) ■ K s {u) . (127) 

Here, the exchange charge g* is defined in terms of the time- component Go(r) of the entan- 
glement vector ( Hill ) by 

g* = J dV G (f) , (128) 

and the sum of both normalisation parameters z a ( 11261 ) - ( 11271 ) yields just the particle 
number (N = 2) 

Zl + z 2 = 2. (129) 

It is true, fixing the value of the normalisation parameters z a ( 11251 ) by the constraints 
( 11261 ) - ( 11271 ) really represents a normalisation condition upon the wave functions V (r), 
namely via the links (1621) - ( 1531 ) between charge densities ^ko(r) and wave functions ip a (r). 
But these conditions still depend on the gauge fields which enter the exchange charge g* 
( 11281 ). Therefore, introducing the "electric" exchange field strength X(f) by 

X(r) = {Xi(r)} = {G 0j (r)} , (130) 

the exchange charge reads in terms of the gauge fields 



9* 



J d 3 (f) B*(r) ■ X(r) - B(r) ■ X*(r) . (131) 



47ras 

Furthermore, the vector field X(f) may be quite generally written in terms of the gauge 
potentials as [l^l 

X{r) = -VSo(r) + iA Q (r)B(r) + iB (r)A 1>2 {f) , (132) 
where the scalar ( Ao ) and the vector ( A1.2 ) potential differences are defined by 

A (r) = — - pU (r) - ^A (f)} (133) 
A h2 {r) = Air) - A 2 (r) . (134) 
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But fortunately, the "electric" exchange potential Bq{t) vanishes identically for the present 
spherically symmetric ansatz (1491) - (152*1) so that the "magnetic" exchange potential B{f) 
becomes proportional to the "electric" exchange field strength X(r): 

X(f) zA (r) • B(r) . (135) 

Using this simplification, together with the functional form ("721) for the exchange potential 
B(f), finally yields for the exchange charge 01281 ) 



3ct s 



j drr 4 A (r)B 2 (r) . (136) 



o 



This result must be inserted into the normalisation conditions ( 11261 ) - ( 11271 ) and thereby 
renders the normalisation procedure an integral part of the eigenvalue problem itself so that 
both parts of the problem must be solved simultaneously. ( For the subsequent numerical 
calculations, an iterative method will be applied where the normalisation ( 11261 ) - ( 11271 ) is 
repeatedly carried out on each iteration step ). For the exact subset of solutions (n\ = ri2 ) of 
section IIX[ the exchange charge g* vanishes because the potential difference A (f) ( 11331 ) 
becomes zero; and thus, the normalisation conditions ( 11261 ) - ( 11271 ) are simplified to z a = 1. 

Finally, it remains to be shown that the postulated normalisation conditions ( 11221 ) are 
consistent with the asymptotic Coulomb form ( 11061 ) of the electrostatic potentials ^A (r). 
For this purpose, rewrite the electrostatic Poisson equations ( "9*151 ) - ( "9*5*1 ) in terms of the 
effective charge densities ^ a %o(r) as 

A (1 U (f) = 47ra s (1) / (r) (137) 
A^A (f) = 4W 2) Z (f) • (138) 

(For the spherically symmetric approximation the effective charge densities ^ a 'lo{r) are here 
to be replaced by their average values over the solid angle Atc). Now, integrating these 
Poisson equations ( 11371 ) - ( 11381 ) over all three - space ( < r < oo ) yields by means of 
GauB' integral theorem ( a = 1, 2 ) 

lim ^1 d 3 fWio(f) _ (139) 

so that the required asymptotic Coulomb form ( 11061 ) of the electrostatic potentials ( a l4o(r) 
is guaranteed just by the normalisation conditions (|122l ) ! (For the topological character of 
the electromagnetic and exchange charges see ref [10( )■ 
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XII. ENERGY FUNCTIONAL 

For the calculation of the energy levels ( Et, say) of the two -electron atoms it is not suf- 
ficient to solve only the mass eigenvalue problem (EH) - ( I1U5I ). because the mass eigenvalues 
M a , a = 1,2, do not immediately determine the energy levels Et- The latter quantity is 
rather fixed by spatial integration of the total energy density ^^Too 

Et — J d 3 f ( T >r o(f) , (140) 

where, in a relativistic theory, the energy density ( T Xr 00 (r) is the time -time component of 
the corresponding energy - momentum density ^T^. This total density is found to consist 
essentially of three parts: the matter density ^T^-, the gauge field density ^T^ v and the 
interaction density due to the external source ^j^: 

^ = ^ + «%, + . (141) 

Accordingly, the total energy Et ( 11401 ) is also built up by three contributions 

Et 

with the self-evident identifications 

E D 
Et 

Ees 

Here, the gauge field density turns out to consist of two subdensities and ^ C ^T^ V 

which are due to the electromagnetic modes F a ^ v (a = l,2)and to the exchange modes 
G^, resp., cf. equation ( II 101 ) : 

(G %, = {K K„ - . (146) 

Naturally, the electromagnetic part appears as the sum of electric and magnetic con- 

tributions ^Tfu, and ^T^, i.e. 

= {c %u + ■ (147) 



Ed + Eg + E e 



d 3 f( D >T 00 (f) 
d 3 f( G >T 00 (f) 
d 3 f ( es >Too(f) . 



(142) 



(143) 
(144) 
(145) 
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A similar splitting does apply also to the exchange density 

{C %u = {h \u + ( %, • (148) 

This splitting of the exchange density ^T^ v is completely analogous to the splitting of its 
electromagnetic counterpart ( 11471 ): namely, the first part ^T^u is the energy - momentum 
density carried by the "electric" field strength X{f)\ and similarly, the second part 
refers to the "magnetic" exchange field strength Y(r). Furthermore, both electromagnetic 
contributions ^T^ v and ^"^T^ are the sums of the mutual - interaction constituents ^T^y 
and ^T^y and of the self- interaction parts and ^ m T^ v : 

(%, = {e % v + (% v (149) 
= + {m % v , (150) 

see ref Since each of these subdensities generates a contribution to the gauge field 

energy via the recipes 0143H - 0145H . the total energy 0142H finally appears as a sum 
of eight contributions 

E T = E D + (E R - E c ) + 4? 

= E D+ (4°> + 4») - E^ - + 4 e) (151) 

= e d + (4) + 4) + 4 m > + 4 m ) - - 4)) + 4 e) • 

The following inspection of all of these energy contributions will yield further insight into 
the logical structure of RST. 

A. Matter Energy Ed 

First, consider the matter energy Ed ( \14'S\ ) which is the corrected sum of both mass- 
energies M a c 2 |3] 

E D = E D{1) + E D{2) = ( Zl ■ M lC 2 - AE m ) + (S3 ■ M 2 c 2 - AE D{2) ) . (152) 

The necessity of correction terms AE-D( a ), a, = 1, 2, is immediately plausible because any mass 
eigenvalue M a (a = l,2)contains the interaction energies of each electron with the other 
one and with the external source; and therefore the simple sum of mass -energy equivalents 
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M a c 2 would embrace not only the sum of rest mass energy and kinetic energy as expected, 
see ref. [10], but also the interaction energies with the other particle and with the external 
source (i.e. nucleus). This is the reason why those mutual interaction energies AEj)( a ) have 
to be subtracted from the individual mass eigenvalues z a ■ M a c 2 in order to obtain the matter 
energy Ed ( 11521 ) exclusively as the sum of the kinetic energies of both particles and nothing 
else. 

Concerning the detailed structure of the subtracted terms Ai? D ( a ), one finds them to be 
a sum of three contributions refering to the external (es) and internal interactions, where 
the latter split up into the electrostatic (e) and exchange (h) types: 

AE D[a) = z a ■ M^c 2 + z a ■ M^c 2 + z a ■ Mf V (153) 
[a = 1,2) . 

The individual contributions are defined by ^| 

c U (f) • {a) k Q (r) (154) 

%(f) ■ <%>(f) (155) 

%(r) ■ <%(i0 (156) 

he f 

- J d 3 f{B (r)h (f + B*(r)h* (f} . (157) 

Clearly, the exchange contribution ( 11571 ) of electric type vanishes for the present para - states 
because here both fields B (r) and h (r) are zero, see the discussion below eq.(l77l). 





M M C 2 = 


-he J d 3 ; 




• M^c 2 = 


-he J d 3 ; 




■ M^c 2 = 


-he J d 3 ; 


Zl 


■ M^c 2 = 


z 2 -Mt ] c- 



B. External Interaction Energy 

A further simplification of the total functional E^ ( 11521 ) occurs in connection with the 
external contributions M^c 2 ( 11541 ): and it is worthwhile to elaborate this point in detail 
because it elucidates the mechanism of electrostatic RST interactions. First, observe here 
that the sum of both external contributions ( 11521 ) brings into play the total charge density 
jo(r): 

2 r 

z a ■ M^c 2 = -he / d 3 f (cx U (f) ■ j (r) (158) 
0=1 J 

jo(f) = {1) Hf + {2 Mr) = ho(f) + 32o(f = - {Mf + j\fj) . (159) 
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But the point with this total charge density jo (r) is now that it generates the total electro- 
static potential A (r), cf. (HBj) 



A (r) = (1 U (r) + (2 Uo(r) 



(160) 



via the ordinary ( i.e. Abelian ) Poisson equation 



d ' 2 d > » - 

cV + ro> Mo(r)=4vra s Jo(r). 



(161) 



Indeed, this claim is easily verified by simply adding up both generalised (i.e. non- 
Abelian) Poisson equations (H23) and (EHl)- Obviously, this interesting effect says that 
an Abelian structure arises for the total objects emerging from the sum of the non - Abelian 
electromagnetic objects! A further demonstration of this effect concerns the electric field 
strengths E a (r) = { {a) F 0j (r}} being originally defined as 

E^f) = -V (1 U (f) - % \B (f)B*(f) - B* (r)B(r) 



E 2 (r) = -V {2) A (r) + 



B (r)B*(r)-B*(r)B(r) 



(162) 
(163) 

Thus, adding up both equations 01621 ) and 01631 ) yields again the ordinary ( i.e. Abelian ) re- 
lationship between the total field strength E(r) and total potential A (r): 



E(r) = £?i(r) + E 2 {r) = -VA {r) . 



Presuming a similar relationship for the external objects 



(164) 



^ex = -VK4 (? 



(165) 



enables one to reexpress the external interaction energy 01581 ) in terms of the electrostatic 
fields strengths E{r) (EM) and E cx ([Ml) 

2 



he 



a=l 



47T«c 



d 3 fE ex -E(r) 



(166) 



This, however, is now a pleasant result because it yields a cancellation of the external 
gauge field energy E es ( 11451 ). Indeed, the external interaction density 

{e %, = ^ | {ex) F, x F„ x + (ex) F, A F/ - ^g^F ffX F° x } (167) 
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is obviously a bilinear construction of the external field strengths ^ ex F^ v ( I11UI ) and the total 
field strength F^ u ( = F 1 ^ + F 2 ^); and this arrangement transcribes to the electric part 
Eet* of the external interaction energy E es ( 11451 ) in the following way 

E£ = ^Jd s rE ex .E(?) (168) 
(i(f) = {^(f)} = {F 0j (f)} , etc) . 

Comparing this to the previous result ( 11661 ) verifies the claimed cancellation 

2 

-Y,Sa-M^c 2 + E^ = (169) 

o=l 

so that the sum of matter energy Ed ( 11521 ) and external interaction energy E^ simplifies 
to 

2 

Ed + J3« = J> " ( M « - M«) c 2 . (170) 

0=1 

Therefore, the final form of the total energy Et ( 11511 ) looks as follows: 
E T = £ *a ■ M a c 2 + (e£ + 4 e) - E K ■ M^A + (e™ + E™) - + . 

0=1 V 0=1 / 

(171) 

This result states that the sum of mass energies is to be corrected by the electrostatic energy 
( first term in brackets ), by the magnetostatic energy ( second term in brackets ), and by the 
exchange energy ( third term in brackets ) . 



C. Internal Electrostatic Interaction 

Naturally, there arises the question whether or not a similar cancellation like those for the 
external contribution in ( 11691 ) can occur also for the internal electrostatic gauge field energy, 
i.e. the first term in brackets on the right-hand side of equation ( 11711 ). The self- energy 
E^ and the mutual- interaction contribution E^ of this term were previously identified as 

4 C) = - g^'M / d¥ + (172) 

4 0) = -^-K p {u) [ d 3 f E^r) ■ E 2 {r) . (173) 
47ra s J 
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Consequently, integrating here by parts does not result in the same cancellation mechanism 
like that for the preceding external case ( 11691 ). Instead, the following relationship between 
the electrostatic gauge field energy E^' and the mass -energy M^c 2 ( 11551 ) - ( 11561 ) holds: 

4 e) = / d 3 f (°>T 00 (r) = 4 C) + 4 e) = i £ z a ■ M^c 2 - N*c 2 . (174) 

J a=l 

The mass - energies M^c 2 read in terms of the wave amplitudes R±, S± 

z x . M[ e) c 2 = -he J dr r 2 {R% + R 2 _ } ( 2 U (r) (175) 

z 2 ■ M^c 2 = -he j dr r 2 {S 2 + + S 2 _ } (1 U (r) , (176) 

and the additional energy content N*c 2 due to the non- linear structure of the non-Abelian 
theory is given by 



AU 2 = ^-^ / <l,VA (1 (r) 



— - Ao(r) 



B 2 (r) . (177) 



Evidently, the electrostatic mass -energy ( 11551 ) - ( 11561 ) of both electrons does not cancel 
against the gauge field energy for the total energy E T ( 11511 ) but, on the contrary, 
has still to be subtracted from the sum of mass eigenvalues together with the additional 
non -linear term N*c 2 . Thus, the total energy E^ ( 11711 ) appears now as 

E T = J2^- M a c 2 -^5>- M^c 2 + nA + (4 m) + 4 m) ) - (4 h) - 4 S) ) . (178) 

a=l \ a=l J 

Clearly, this is again a plausible result because the sum of mass eigenvalues ( first term on the 
right - hand side ) counts twice the interelectronic interaction energy of the electrostatic type 
which therefore has to be subtracted once (second term on the right-hand side). Observe 
also, that this electrostatic interaction energy, being defined by the mean value of expressions 
( 11551 ) - ( 11561 ). contains no manifest self - coupling of the form ■ ^'Aq. But an indirect 
self - coupling does clearly exist because the a-th potential w A${r) is partly generated by 
the 6-th wave function ipj, (a ^ b), see the electrostatic Poisson equations 



D. Internal Magnet ost at ic Interaction 

After the role of the electrostatic interactions has been clarified, one may turn attention 
to the somewhat different pattern of their magnetostatic counterparts. The mutual and 
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self - interaction parts of this type of interactions were identified in terms of the magnetic 
fields H a {r) as Irj 

4 m) = -^K s (u) [ d 3 f {Hl{f) + # 2 2 (f)} (179) 
4 m) = "t^^pH / d 3 f^(f) • H 2 (r) . (180) 

47T« S J 

It is true, this looks quite similar to the electrostatic counterpart 01721 ) - 01731 ) but observe 
that for the magnetostatic interaction energy, in contrast to the electrostatic case, there is 
no need of removing it from the matter energy Ed ( 11521 ). Another difference between both 
types of interactions is the non-Abelian construction of the magnetostatic fields H a (r) 

H a (r) = {^W(r)} = jl^V a) iV} (181) 

by means of the potentials A a (r), B(r); see ref 

H^r) = V x Ai(r) - iB(f) x B*(f) (182) 
H 2 (r) = V x A 2 (r) + iB(r) x B*(f) . (183) 

This differs from the analogous situation with the electrostatic gradient fields E a {r) ( 11621 ) - 

(EEB1) 

E a (r) = -V H4 (r) (184) 
'a = 1,2 N 



where the non- linear terms become zero on account of the vanishing "electric" exchange 
potential Bq{t)\ 

Clearly, this fact renders more complicated the mutual and self- interaction energies Ej^ 
and Er of both particles 

4 m) = ~^ K M J d'mir) ■ H 2 (r) (185) 

4 m) = -^-K p (u) [ d¥ \Ht(r) + #|(f)| . (186) 

But on principle, the procedure is here the same as for the electrostatic case, i.e. one inserts 
the functional form of the magnetic fields H a (f) (due to the isotropic ansatz) 

H x {r) = 2costf [Ai(r) - r 2 B 2 (r)] ■ e r - (^(r)) • (187) 

— * sin d 

H 2 (r) = 2 cost? [A 2 (r) + r 2 B 2 {r)] ■ e r - —j- {r 2 A 2 {r)) ■ e# , (188) 
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and then one finds the magnetostatic energy functional ( 11851 ) and ( 11861 ) by partially in- 
tegrating and use of the magnetostatic Poisson equations ( 11001 ) and ( 11011 ) to be of the 
following form: 

pH 2 



+ 2r 2 5 2 } 



M m) = -hcK s e~ 2u 
r 3 



2r 2 B 2 



J drr^A' 2 -R + R_-A[-S + S_} + ^-^ J dr r 4 B 2 [ A 1)2 (l - r 2 A li2 ) 

(189) 

J drr 3 {A[-R + R„-A' 2 .S + S-}--£p^ J dr r 4 5 2 { A 1)2 (l - r 2 A lj2 ) 

(190) 



Despite the formal similarities of these magnetostatic energy contributions E^\ E { ™ ] (MB) 
and ( 11861 ) with their electromagnetic counterparts E^ , ( 11721 ) and ( 11731 ) there is an 
important difference which refers to the occurrence of direct self - interactions for the mag- 
netostatic case; see the remarks below equation ( 11781 ) for the corresponding electrostatic 
situation. The point with the magnetic case refers here to the occurrence of modified mag- 
netostatic potentials A' a (r) which are linked to the original A a (r) (1701)- (I7H) in the following 
way: 



A[(r) = K s (u) ■ A 2 (r) - K p (u) ■ A^r) 
A' 2 (r) = K s (u) ■ A 1 (r) - K p (u) ■ A 2 (r) 



(191) 
(192) 



This link clearly reveals the existence of direct self- interactions of the magnetic type since, 
e.g., the second modified potential A' 2 in ( 11891 ). on the one hand, combines with the wave 
amplitudes R± of the first particle, but on the other hand, it contains the first potential 
A\[r) which itself is generated just by the first particle (a = l)itself. In this way, the first 
particle acts back to itself: 



R± => A 1 (r) => A' 2 (r) 



A' 2 ■ R+R_ 



E 



(m) 
R 



(193) 



This reasoning holds analogously for the first magnetostatic mode A[. But clearly, an even 
more manifest self - interaction is established by the proper self-energy E { ™ ] (HSU) where 
the first (second) magnetostatic potential A\ (A' 2 ) directly combines with the first (second) 
current intensity k\ (/c 2 ), cf. eq.s ( 11031 )- ( fT04l ). 
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E. Exchange Interactions 



The last contribution to the total energy is the exchange energy Ec, which is the energy 
content of the exchange modes G^ u as specified by ( 11481 ). i.e. 

E c = [ d 3 f (°>To (f) . (194) 



Since the energy - momentum density ( 11481 ) of the exchange modes splits up into an 

"electric" part and a magnetic part ^^T^u, this splitting transcribes to the exchange 

energy E c 



E c = 4 h) + 4 g) (195) 
4 h) = J d 3 f ^>Too(r) (196) 

= J d 3 f^% (r) . (197) 

Concerning the "electric" part E c h \ it is true that both the exchange density ho(r) and the 
associated scalar exchange potential Bq{t) do vanish identically for the para -helium states 
x Sq. But nevertheless, the exchange field strength X(r) is non- trivial, cf. ( 11351 ). Therefore, 
the corresponding exchange energy EjP of the "electricity pe is non -zero and reads, quite 
generally, in terms of the exchange field strength X{f) [13j 

,(h) _ he f ^3-.^* 
47ras e_2n 



#c = ~ A =^7 / d"fX*(r) ■ X(r) , (198) 



i.e. by use of the exchange field strength X(f) (|T 

4 h) = |^I / drr*Al(r)B\r). (199) 

The final contribution to the total energy Et ( 11781 ) is the exchange energy of "magnetic" 
type H 

E c ] = A „ h °^2u I d 3 rY*(r) ■ Y(r) (200) 



47ro;se 



Decomposing here the "magnetic" exchange field strength Y(r) with respect to spherical 
polar coordinates as usual 

Y(r) = Y r e r + Y#e# + Y v e v , (201) 
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with components 



Y r = -2 sin i? e ttp B(r) 



1 - ir 2 A li2 (r) 



Y« = -costfe** 



l_d_ 

r dr 



[r 2 £?( 



?C 



1 d 



r dr 



(202) 

(203) 
(204) 



ultimately yields by means of integrating by parts and use of the exchange Poisson equation 



i una i 



E, 



(g) 



-he 



[ drr 3 B(r) [R + ■ + R_ ■ S+] l — [ drr A B 2 (r) 

J "se lu J 



A 2 (r) + A li2 (r) 



2r 2 B 2 (r) 



(205) 

With this result, the total energy Et is completely fixed in terms of the wave amplitudes 
R±(r), S±(r) and the gauge potentials ^A (r), A a (r) and B(r). Its value upon the solutions 
of the energy eigenvalue problem, which consists of the mass eigenvalue equations and the 
Poisson equations, will yield the desired energy levels of the ground state E^' 1 ^ and the first 

(1 2) 

excited state E™ ' in the next section. 



F. Special Case: ns 2 

As a check for the more general results concerning the states n\s n 2 s ^o, return to 
the subset of exactly isotropic states ns 2 ^ for which the total energy E T was previously 
determined as jl^ : 

E T = 2M"c 2 - E^ 1 + 3^ m) , (206) 
with the electrostatic energy E^ being given by 



E® = -He J drr 2 [R 2 + (r) + R 2 _(r)] 



o(r) , (207) 



and its magnetic counterpart E^ by 



#H = ^ c J J & rr ZB(r)R + {r)R^{r) - - *_ 2u J drr 4 B 3 (r) [l - r 2 B{r)] \ . (208) 



o o 
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Observe the prefactor of three for the magnetostatic contribution ( 12061 ) to the total 
energy . It was shown in the preceeding articles 0, Q] that this threefold weight of E^ 
relative to E^ is due to the exact isotropy of the states ns 2 ^So, where each of the three 
coordinate axis contributes the same amount to the spin -spin interaction energy. 

The special form ( 12061 ) of the energy functional Et is immediately obtainable from 
the more general result ( 11781 ) by use of the former identifications ([781) - (18*31) . They put 
the exchange charge g* ( 11311 ) to zero because the "electric" exchange field strength X(r) 
vanishes as does the electrostatic potential difference A (r) ( 11331 ). As a consequence, both 
normalisation parameters z a ( 11261 ) - ( 11271 ) adopt unity value, z a =>- 1. In turn, this result 
reduces the general sum of mass eigenvalues occuring in the energy functional ( 11781 ) to 
the special form ( 12061 ): 

2 

z a ■ M a c 2 2M"c 2 . (209) 

a=l 

Next, pay attention to the vanishing of the non- linear term N*c 2 ( 11771 ) which is due to 
the vanishing electrostatic potential difference A (r) ( 11331 ): 

N*c 2 . (210) 

Thus, the remainder of the electrostatic interactions equals half the sum of mass equivalents 
Ma^ ( 11551 ) - ( 11561 ) as being specified by the results ( 11751 ) - ( 11761 ). By the identifications 
(ESI) - (ESI) the latter two equations adopt the following special form 

1 2 r 

=> -he / drr 2ip) Ao(r){R 2 + + R 2 _} . (211) 

This is just the required result for the electromagnetic interaction energy E ( l ] flUTI) due to 

the states ns 2 x Sq. 

Next, consider the general mag netostatic contributions E^ ] ( 11831 ) and E^ (JHM). 
Again, applying to these expressions the identifications ( 1781 ) - (EH]) simplifies them to 

£H -ihcK v (u)e- 2u jy drr 3 R + (r)R4r)B(r) - - *_ 2u J drr 4 B 3 (r) [l - r 2 B(r)] j 

(212) 

j?H => ~^hcK s {u)e- 2u | j &rr 3 R + (r)R_(r)B(r) - — ^ j drr 4 B 3 (r) [l - r 2 B(r)] J . 

(213) 
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Therefore, for the exactly isotropic states ns 2 ^So, one adds up both equations ( 12121 ) and 
( 12131 ) and thus finds the required special form ( 12081 ). 

To conclude matters, call your attention to Eq ( 11991 ) and E^ ( 12051 ). the yet uncon- 
sidered exchange energy contributions of "electric" and "magnetic" type to the total energy 
-&r H178I ). Clearly, the "electric" contribution Eq vanishes because of the vanishing elec- 
trostatic potential difference A (r) ( 11331 ): 

E^ . (214) 

Up to sign, the magnetic contribution E^ ( 12051 ) adopts the double value of the magneto- 
static contribution ( 12081 ): 

4 g) -2E^ ] . (215) 

So, in summary, the special form ( 12061 ) of the energy functional Et for the exact isotropic 
states does indeed emerge from the general shape ( 11781 ). namely by substituting therein the 
special results ( 12091 ) - ( 12151 ) referring to the states ns 2 1 Sq. 



XIII. NUMERICAL RESULTS 



For a comparison of the present RST results to the experimental data and to the cor- 
responding predictions of other theoretical approaches, observe that the literature mostly 
presents ionisation energies ( jjh , say) which refer to the states of the helium- like ions, 
where one of both electrons occupies the groundstate (Is) and the other one is in an excited 
state (ns). 

Hence, comparing the various theoretical predictions for the chosen energy difference 



Hheo 



£i_2 4- E theo \ ls 2s 1So - -E t heo|i s 2 i So (216) 



among each other and all of them to the experimental values A exp .Ei_2, has to be preceded 
by transcribing the ionisation energies J^l 2 ) ^° ^ ne desired energy difference ( 12161 ). This 
can be easily done by simply observing that the concept of "ionisation energy" refers to the 
energy difference between a final state and an initial state, where in the final state one of 
both electrons is free, having energy Mc 2 , and the other one occupies a bound state with 
energy -E t heo(i)- On the other hand, in the initial state both electrons do occupy a two- 
particle state of energy i\heo(2) - Thus, the theoretical ionisation energy Jy[eo(2) ^ s & ven by 
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J$&) = + <L (1) - ■ (217) 

Putting now n = 1 and n = 2 , and subtracting the corresponding equations results in 
the desired energy difference A thco £ , 1 _ 2 : 

Athco^l-2 = ^ t io(2) ~ J th£(2) " ( 218 ) 

Here the two - particle ionisation energies J t ^L ^ and J^eo(2) ma y t" e taken from the vari- 
ous theoretical approaches in the literature |llL Il2|. w hereas for the one -particle energies 
-^thL(i) (ETZI) the generally accepted values of ref. [l^f are suitable. 

It must be stressed that we restrict ourselves to a point - like and fixed nucleus ( of charge 
number 2 < z cx < 100 ) so that the effects due to mass polarisation and nuclear size are 
neglected. Therefore, only the other contributions to the Lamb shifts of the one -and two- 
electron configurations are taken into account. This can be made more concrete by represent- 
ing the energy difference A t h eo -&L-2 ( 12181 ) as the sum of its "semiclassical" part Atheo-^1-2 
and the self- energy correction 5 s Ei- 2 

Atheo-E'i-2 = A theo _E , 1 _ 2 + o~ s Ei-2 ■ (219) 

Here, the semiclassical contribution A t h eo -Ei-2 is given in terms of the semi - classical ioni- 
sation energies ^30(2) as 

A thcoJ Ei„2 = J t teo(2) ~ ^theo(2) ' ( 220 ) 

where the ionization energies themselves are linked to the one- and two-particle energy levels 

7W Q nJ p(M 



^theo(i) and E Lo(2) through 



f(l.n) 
t7 theo(2) 



Mc 2 + 



E^\ 2) ■ (221) 



The energy eigenvalues M^c 2 of the one -particle Dirac - Coulomb problem can be looked 
up in any textbook about relativistic quantum mechanics, e.g. j^j]: 

£(n) = M (n) 2 _ M ° 2 (222) 

^theo(i) - M * c - g | (2cxQs)2 • 



(n-l+x/l-(«exO;s) 2 ) 

Accordingly, the self- energy correction 5 S E 1 _ 2 in the standard approaches is composed of 
the QED corrections ( Lamb shift ) of the ionisation energies: 

8 S E^ 2 = 8 S J t H i2) ~ <^theo(2) • ( 223 ) 
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The conventional Lamb shifts ^^00(2) °f the ionisation energies, with and without mass 



polarisation and nuclear size effects, are tabulated in ref. [11], so that a detailed comparison 
with the corresponding RST predictions becomes now feasible. 

A. Semiclassical Results 

For the semiclassical approach, one neglects in the conventional theories the Lamb shift 
which in RST corresponds to putting the self- interaction parameter u to zero. Under 
these presumptions, Table I presents a comparison of the theoretical predictions for the 



semiclassical energy difference AE\^2 due to the relativistic 1/Z-expansion method |U| and 
RST. Here the RST procedure consists in solving numerically the mass eigenvalue problem 
of Sect. |3 and taking the value of the energy functional E? ()178|) upon the numerical 
solutions. 

The experimental values A oxp i?i_ 2 jl^ are also displayed in Table I, though it is clear 
that their coincidence with the theoretical predictions can be achieved only after including 
the self-energy effects (see below). In order to facilitate the analysis of the relationship 
between theory and experiment, the numerical content of table I is illustrated by fig.l 
through displaying the relative deviations (A) from the 1/Z expansion method: 

o o 

Arst E1-2 —Ayz E1-2 



iRST 



A V z E 



(224a) 



1-2 



A exp = KxpEl - 2 - Al/Z El - 2 . (224b) 
Ai/z E1-2 

The semiclassical RST predictions Arst E1-2 (7th column of Table I) deviate from 
the experimental values A cxp _E'i_ 2 by 1,4% for z ex = 10 (neon) and decrease down to 
0,2% for z ex = 42 (molybdenum). This effect of improvement of the RST predictions for 
increasing nuclear charge number z ex is observed also in other situations, see ref. 131 ]; but in 
view of the more precise 1/Z-predictions (4th column of Table I) this result requires some 

o — — — — 

explanation. The point here is that the energy difference Arst E1-2 (jZZDj) is composed 

0(1,1) 0(1,2) 

of the two ionization energies Jrst(2) an d </rst(2) which therefore both do influence 

0(1,1) 

the numerical value for Arst Ei-2- Here, the ground-state ionization energy J RST ( 2 ) 
is calculated exactly and is found to be in sufficient agreement with the corresponding 



=(i,i) 



=(1.2) 



1/Z-predictions J 1/2(2)1 see re f-|l2|- However, the ionization energy Jr ST ( 2 ) of the first 
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excited state Is 2s 1 S'o is calculated by means of the spherically symmetric approximation 
(see the remarks above the eigenvalue equations ()91 )) -(|94 )) ): and this approximation 
introduces an additional (albeit artificial) deviation which mainly is responsible for the 
relatively large deviations of the semiclassical RST predictions from the experimental values. 



Table I 

This effect may be exemplified numerically, e.g., for nuclear charge number z ex = 30 
(zinc), cf. Table I. Here the 1/Z prediction amounts to Au z E\-2= 8 955.4 eV and thus 
deviates from the experimental value A exp i?i_2 = 8 950.2 eV by 0.06%. In contrast to this, 

o 

the corresponding RST prediction is A RST Ei-2= 8 912.2 eV and therefore has a deviation 
of 0.4%. However, this relatively large deviation is almost completely due to the ionization 

o(l,2) 

energy Jrst(2) = 2 957.5 eV of the excited state, which deviates from the corresponding 

o(l,2) 

1/Z value J i/^(2)— 2 915.67 eV by 1.4%. On the other hand, the RST prediction for the 
ionization energy (|221j) of the groundstate is Jrst(2) = H 869.7 eV and therefore deviates 

o(l,l) 

from the corresponding 1/Z value J \/z(p,)— H 871.03 eV by only 0.01%. Consequently, it is 

°(1>2) 

just the relatively inaccurate RST prediction J R g T ( 2 ) for the excited state (being caused by 
the spherically symmetric approximation) which is responsible for the large RST deviation 
from the experimental value. This result concretizes the necessity of developping some 
methods to solve exactly the two-dimensional RST eigenvalue problems in order to avoid 
the use of the spherically symmetric approximation. Alternatively, when one wishes to stick 
to the spherically symmetric approximation, one might also look for a better class of trial 
functions than those given by equations (}4*9"j) - (p)2"j) . 

fig.l 

B. Self-Interactions 

The inclusion of the electronic self-energy works somewhat different in RST as in the 
conventional theory: The RST self-interactions constitute an integral part of the theory, 
namely in specifying the fibre metric K a p (190 jl which must be fixed in any way, no mat- 
ter whether one wants to include or neglect the self-interactions. In contrast to this, the 
self-interactions of the conventional theory are added by hand to the semiclassical results 
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in form of perturbative corrections due to quantum electrodynamics. The latter theory is a 
rather different framework which however is not able to yield immediately the desired energy 
levels of atoms and molecules, so that one is forced to resort to a two-track procedure (i.e. 
calculating semiclassical results and adding QED corrections). A typical example for this 
is the relativistic 1/Z-expansion method of ref. |h| whose QED corrections have been used 
also for the all-order technique in relativistic many-body perturbation theory 12 1. 

Thus, it becomes interesting to oppose these rather different approaches to each other 
by comparing the QED-corrected semiclassical predictions Ayz -&1-2 of ref.s [ill Q to 
the analogous RST results Arst E1-2 with non-zero self-interaction parameter u. Here the 
precise value of u is adopted as u = 0.03052, because for this value of u the electronic 
ground-state interaction energies and ionization energies have been found to coincide ex- 
cellently with the experimental facts for larg ref.s jiol Q]. Clearly one wishes 



to know whether the self-interaction parameter u must adopt a universal value, or whether 
perhaps it depends upon the quantum state to be considered. Table II and fig. 2 present 
a comparison of the RST and 1/Z predictions Arst-Ei-2 arid Ai/z-Ei-9 with inclusion of 



ft 



the self-energy effects. The 1/Z data are directly taken over from ref. [Uj with neglection 
of mass polarization and nuclear size, and the corresponding RST data are obtained by 
numerically solving the full RST mass-eigenvalue system (l9*Tl)-( [TUo1 ) with u = 0.03052. 

Concerning the QED corrected 1/Z predictions A 1 /^i? 1 _ 2 ([11]) (fourth column of Ta- 
ble II), the inclusion of the QED self-interactions yields an acceptable agreement with the 
experimental data up to medium nuclear charge numbers (z ex < 50), see also fig. 2. Here 
the inclusion of the QED self-energy actually generates a considerable improvement. In 
contrast to this, the inclusion of the RST self-interactions (by taking u = 0.03052) results in 
an ambiguous shift of the energy difference Arst-E'i-2 (( 12181 )). see seventh column of Table 
II. For lower z ex (< 28), the inclusion of the RST self-interactions shifts the RST predic- 
tions Arst-E'i-2 (( 12181 )) closer to the experimental data, whereas for the higher z ex (> 28) 
the corresponding RST predictions move apart from the experimental values. However, the 
semiclassical tendency of better appoximation for the higher nuclear charge numbers does 
persist also after inclusion of the self-interactions. 

fig.2 
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In this sense, RST seems to represent some kind of high-energy approximation to the real 
world. Indeed, this feature of RST is even more manifest for the ground-state ionization 
energies of the helium-like ions, for which the RST predictions are more precise than the 
predictions of the other theoretical approaches in the literature [l^ . Clearly, such a result 
hints at the unexhausted potentiality of RST, where it seems possible to further improve 
the RST predictions by finding approximations to the exact solutions which supersede the 
present trial functions in combination with the variational technique. 

Table II 
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all data in eV 



Table I: Semiclassical Energy Differences A Ei-2 

The semiclassical results of the l jZ expansion method ^ | (fourth column) for the excitation 

o o 

energy A E1-2 are compared to the present RST predictions Arst E1-2 according to equation 
(|22U|) . The relatively large deviations of the RST results (seventh column) from the experimental 
values (last column) are attributed to the anisotropy effect, which is neglecte d by the spherically 
symmetric approximation. The experimental values are taken from reference 
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all data in eV 



Table II: Energy Differences AE1-2 with Self- Interaction 

The semiclassical predictions for the energy difference A£i_2 become improved considerably 
for the conventional 1/Z expansion method but only partly for RST by including the self-energy 
effects. The moderate precision of the RST predictions is a consequence of the use of (i) the 
spherically symmetric approximation and (ii) the simple one-particle trial functions (EH) - (EH) • 
It is only the use of these very rough approximations which devaluates the RST predictions down 
to the level of the simple Hartree-Fock approach. 
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Fig.l: Relative Deviations A (|^4 )- |M^ 

The relative semiclassical deviations from the exp erimental values (t) are larger for the RST 
approach (♦) than for the 1/Z expansion method because the spherically symmetric 
approximation of RST induces an extra (albeit artificial) deviation in addition to the missing of 
the correlation energy beyond Haxtree-Fock. For high nuclear charge numbers z ex , the RST 



predictions (♦) approximate the 1/Z predictions (horizonta 



I •:, r >> i ! coincidence of RST and 1/Z predictions see i;ef. 1^ | 



axis). For this asymptotic 
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Fig. 2: Relative Deviations (\224a\ )- ^224^ ) with Self- Interactions 

The inclusion of the QED corrections yields a considerable improvement of the 1/Z 
predictions (t), see fig.l; whereas the modification of the RST predictions (♦) is ambiguous: 
improvement for low z ex {< 28) but deterioration for the higher z ex (> 28). This is interpreted as 
an artifical effect due to the use of the spherically symmetric approximation in combination with 
the simple one - particle trial functions ( 1491 ) - (E21) . 
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